Bcc 4 He as a Coherent Quantum Solid 
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In this work we investigate implications of the quantum nature of bcc 4 He. We show that it 
is a unique solid phase with both a lattice structure and an Off-Diagonal Long Range Order of 
coherently oscillating local electric dipole moments. These dipoles arise from the local motion of 
the atoms in the crystal potential well, and oscillate in synchrony to reduce the dipolar interaction 
energy. The dipolar ground-state is therefore found to be a coherent state with a well defined global 
phase and a three-component complex order parameter. The condensation energy of the dipoles in 
the bcc phase stabilizes it over the hep phase at finite temperatures. We further show that there 
can be fermionic excitations of this ground-state and predict that they form an optical-like branch 
in the (110) direction. A comparison with 'super-solid' models is also discussed. 



PACS: 67.80-s,67.80.Cx,67.80.Mg 
I. INTRODUCTION 

The bcc phase of 4 He has a pronounced quantum na- 
ture due to the relatively open structure of the lattice. 
Quantum effects are manifested in strong anharmonicity 
of some phonon modes and in the large zero-point ki- 
netic energy of the atomsEl. It is this large kinetic en- 
ergy which is thought to help stabilize the bcc phase 
over the hep phase. In this paper we highlight this na- 
ture of the bcc phase by proposing a new physical model 
for the local atomic motion. For the sake of clarity we 
have reproduced here gome of the arguments and calcula- 
tions already given ina. We propose that in the bcc 4 He 
phase the local excitations of the atoms in their poten- 
tial wells, result in oscillating local electric dipoles. The 
ground-state of these dipoles has the dipoles oscillating 
in synchrony, thereby reducing the dipolar interaction 
energy between them. Solving a mean-field Hamiltonian 
describing these dipoles we find that Bosonic phase fluc- 
tuations in the (110) direction reproduce the spectrum 
of the Ti(110) phonon. 

In the following we further explore the nature of the 
coherent ground-state of the local-modes in the bcc 4 He. 
We show that the bcc 4 He is a unique phase having both 
Diagonal Long Range Order (DLRO) of the solid lat- 
tice and Off-Diagonal Long Range Order (ODLRO) of 
the local dipoles. There is therefore a complex three- 
component order parameter which describes the coher- 
ently oscillating dipoles in each of the three orthogonal 
directions in the lattice. In the ground-state the local 
dipoles form a Bose-Einstein condensate in the zero mo- 
mentum state, and we are able to estimate the ground- 
state energy reduction due to this condensation. This 
estimate compares favourably with experimental results 
and consequently we claim that this condensation energy 
stabilizes the bcc phase over the hep phase. In the hep 
phase we expect no coherence or condensation due to the 
highly isotropic lattice and the geometric frustration of 



the hexagonal symmetry. We also comment about the 
relation of this work to previous work about the 'super- 
solid' concept in quantum solids. 

Additionally, we predict a high-energy optical-like 
mode which has fermionic statistics. This excitation is 
confined to the (110) direction and involves a local 'flip- 
ing' of a dipole with respect to the ground-state. This 
makes the dipole become anti-symmetric (a tt phase dif- 
ference) with respect to the global phase of the complex 
order parameter and aquire Fermi-Dirac statistics. We 
also give analytic expressions for the scattering intensity 
of both the Bose and Fermi excitations along the (110) 
direction. These predictions remain to be compared with 
future experimental data. 



II. GROUND-STATE COHERENCE AND BOSE 
EXCITATIONS. 

The usual treatment of the ground-state and energy 
of bcc He employs variational wavefunctions that aim 
to account for the short-range correlations between the 
atomsa. These correlations arise mainly due to the hard- 
core repulsion between the atoms. The atoms have a 
high zero-point kinetic energy which is given quite accu- 
rately by treating them as independent particles held in 
place by the potential of the neighboring (static) atoms. 
This type of calculation is the " particle-in-cell" approx- 
imation which gives surprisingly good agreement with 
measured thermodynamic properties of the solid phaseS. 
We want to focus here on the effects of the local motion 
of the atoms inside this potential-well on the nature of 
the ground-state. In this approach we would like to iso- 
late the lowest energy excited state of the atom inside 
its potential well, and treat it as a local excitation of the 
lattice. This local excited state consists of a local oscilla- 
tory motion of the atom along a particular direction and 
produces an oscillating electric dipole, similarly to that 
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of the usual Van-der Waals interaction. However unlike 
the case of the Van-der Waals interaction, in which the 
dipolar fluctuations are random, we show that in the bcc 
solid there are local dipoles which are correlated and a 
new ground-state of lower energy is created. 

The potential well of an atom in the bcc lattice due to 
the standard helium pair-potential v(rp, provided one 
can take the other atoms as stationary, can be maped 
along any direction in the lattice. We findcl that in the 
directions normal to the unit cube's faces (i.e. (100), (010) 
etc.) the confining potential well is very wide with a pro- 
nounced double- minimum structure (Fig.l). Solving the 
one-dimensional Schrodinger's equation for a 4 Hc atom 
in this potential, we get a first excited level with energy 
10K, and a wavefunction describing atomic motion with 
an amplitude of ~ lA (in the (100) direction (Fig.l)). 
Based on the above calculation, we shall assume that the 
atoms have a local-mode that is highly directional along 
one of the directions equivalent to (100). Local atomic 
motion along the other directions is assumed to be severly 
restricted due to the higher excitation energies (Fig.l). 
Experimental evidence for the existence of such a "local 
mode" comes from NMRj-measurements which find an 
activation energy of 7±1KqIj, and we propose to identify 
this local-mode with the highly directional motion of the 
atoms in the normal directions. 

Using this identification we can now estimate the size 
of the local electric dipole moment that can be created 
by this local and highly directional atomic motion. As 
the atom moves this instantaneous local electric-dipole 
is created due to the electronic cloud and the ion be- 
ing slightly displaced relative to each other. The electric 
dipole moment due to mixing of the lowest \s) and \p) 
electronic-levels of the 4 He atom, is given from perturba- 
tion theory as 

V> = |s) + A \p) ^E Q ~(^ \E\ $ - (a \E\ s) ~ A 2 (p \E\ p) 
=> A 2 ~ 7/2.46 • 10 4 ~ 0.00284, A ~ 0.0168 (1) 

where |s) and \p) stand for the ground-state and first 
excited-state of the 4 He atom, A is the mixing coefficient 
and (p \E\ p) ~ 2.46 • 10 4 K is the excitation energy of the 
first atomic excited-statecl. This small estimated mixing 
gives the magnitude of the induced dipole moment as 



instantaneous dipolar interaction energy for each of the 
three orthogonal directions, is given by 



= e{ip \x\ ip) ~ 2eA (s \x\p) ~ e • 0.03A 



(2) 



where (s \x\ p) ~ 0.9A . The estimation of the mixing A 
and the dipole-moment serves to set an upper bound 
on the magnitude of this effect, since we assumed that 
the entire excitation energy Eq is converted to a local 
electric dipole. 

It is possible to show that the lowest energy of a corre- 
lated dipolar array in the bcc lattice preserves the sym- 
metry of the bcc unit cell along one of the symmetry axes. 
In such a case it can be easily shown that there will be 
no contribution to the dipolar interaction energy from os- 
cillating dipole moments which are orthogonal, and the 
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3 cos 2 (fi ■ (r - ri)) - 1 



(3) 



where the sum is over all the atoms in the lattice, Vi 
being the instantaneous coordinate of the i-th atom. For 
oscillating dipoles with random phases, the average in- 
stantaneous interaction energy (|3|) summed over the lat- 
tice would be zero. However, the energy of the dipoles 
can be lowered by correlating the phases of the oscillating 
atoms. Since the direction of the local dipole shows the 
instantaneous direction of the motion or displacement, a 
state where all the dipoles point in the same direction is 
just a uniform motion or translation of the entire lattice. 
We therefore have to look for symmetric arrangements 
with respect to the number of up/down dipoles, such as 
is shown in Fig. 2. This is the lowest energy 'antifer- 
roelectric' configuration with the periodicity of the bcc 
unit cell. We have shown this arrangement for individual 
dipoles oriented along the (001) direction, but they are 
similarly arranged for dipoles along the two other orthog- 
onal axes. The sum in (|3|) for such a configuration with 
a unit dipole is given in Fig. 2. Thus, the ground state in 
our picture has the atoms executing this correlated local 
oscillation along the three orthogonal directions. 

We therefore have, in addition to the usual 
(isotropic) Van-der Waals interaction, highly directional 
(anisotropic) electric dipoles that become correlated so 
that they oscillate in synchrony. This is a state of quan- 
tum resonance where the system oscillates between two 
equivalent up/down arrangements of the ground-state of 
the dipoles (Fig. 2). The total interaction between the 
atoms is now given as the usual second-order (cx 1/r 6 ) 
Van-der Waals contribution that is the result of local- 
dipoles which have random relative phases, and an ad- 
ditional long-range (first-order, oc 1/r 3 ) dipolar interac- 
tion from the correlated part. Dipolar interactions that 
decay as 1/r 3 occur for perfectly correlated oscillating 
dipoles, such as a single electric dipole and its image in 
an adjacent conducting plate. The coherently oscillating 
nearest-neighbor dipoles therefore behave as perfect im- 
ages of each other (Fig. 2), and oscillate with the same 
global phase. 

The correlated oscillating dipoles do not have an av- 
erage static dipole moment, so this is not the case of an 
antifcrroelectric structural phase transition!! The array 
shown in Fig. 2 is simultaneously arranged along the other 
two orthogonal axes. Along each direction the ground- 
state is given as a coherent-state of these local dipoles, 
i.e. has a well-defined phase and an ill-defined occupation 
number. 

We shall treat the dynamics of the correlated dipo- 
lar array as independent of the other degrees of free- 
dom of the lattice. This assumption needs justification 
since there can be phonon modes that will modulate 
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the atomic motion, thereby coupling with the oscillat- 
ing dipolar array. The oscillatory atomic motion induced 
by the phonons will modulate the relative phases of the 
dipoles. Let us look at the ground state of the dipoles, 
taking for example dipolar oscillations oriented along the 
(001) direction (Fig. 2). We now need to consider only 
phonons which will modulate the local motion responsi- 
ble for the oscillating dipoles in this direction. In the bec 
structure, only 3 phonons fulfill this condition: L(001), 
T(100) and Ti(110). Let us calculate the energy of the 
dipolar array when modulated along these 3 directions. 
For a modulation along some, direction k , the dipolar 
interaction energy is given byu: 



X(k) =-| M | 



3 cos 2 (/i • (r - fj)) - 1 



exp [2nik ■ (r 



|ro- 



(4) 



At fc = the interaction matrix X(k) is just the dipolar 
energy (§). 

In Fig. 3 we plot the value of X(k), the energy of the 
dipolar array modulated by the relevant phonons: L(001), 
T(100), and Ti(110), for dipole moment = 1. We 
see that for a modulation by L(001) and T(100) the pe- 
riodicity of X(k) is over a full unit-cell, that is twice 
the periodicity of these phonons. Since symmetric func- 
tions of periodicities Tr/a and 2n/a are orthogonal, so 
are the eigenfunctions of these particular phonons and 
dipole-excitations. The dipole array cannot therefore be 
excited by any of these two phonon. For the modulation 
produced by the Ti(110) mode, the periodicity of X(k) 
is the same as that of the Ti(110) phonon, which can 
therefore couple to the dipole array. We conlude there- 
fore that the coupling of the local modes to the lattice 
excitations is limited to a single phonon mode, justifying 
our assumption that the local modes can be treated sep- 
arately to a good approximation. We shall now calculate 
the dispersion relation of such an excitation by a mean- 
field solution of an effective Hamiltonian. It turns out 
that the only phonon mode of the bec lattice that can 
couple with the dipolar array is in fact the natural exci- 
tation of the dipolar array in the (110) direction. Thus, 
the only elementary (Bose) excitations of the dipole ar- 
ray would be in the (110) direction, in the form of the 
Ti(110) phonon. The description of this phonon is there- 
fore taken into account by our treatment of the dynamics 
of the dipolar array, and will appear as a solution of the 
mean- field treatment. This means that our assumption 
of an effective decoupling between the dipolar and other 
degrees of freedom is justified. 

The Hamiltonian treatment of interacting local excita- 
tions was developed originally by Hopfielda for the prob- 
lem of excitons in a dielectric material. The local excita- 
tions are treated as bosons using the standard Holstein- 
Primakof procedure, and-jthe effective Hamiltonian de- 
scribing their behavior ialj 



H loc = ^(£ + X{k)) (bjb k + X - 

k > 



(5) 



where b^ ,bk are Bose creation/anihilation operators of 
the local mode, X(k) is given in (||) and E is the energy 
of exciting a local dipole out of the correlated ground- 
state. 

The Hamiltonian H\ oc (||) which describes the effec- 
tive interaction between localized modes can be diagonal- 
ized using the Bogoliubov transformation = u{k)bk + 
v(k)tf The two functions u(k) and v (fc) are given by: 



u 2 (k) 



E +X(k) 
E{k) 



+ l),v 2 (k) 



E +X(k) 
E(k) 



(6) 



The result of solving by mean-field thejjeffective Hamil- 
tonian for the correlated dipolar arrayo, is a coherent 
ground-state given byO 



and the energy spectrum is 



E(k) = y/Eo {Eo + 2X(k)) 



(7) 



(8) 



In Fig. 2 we see that the ground-state arrangement has 
the dipoles arranged in alternating planes in the (110) 
direction. As we have shown the only naturally occur- 
ing Bose excitations of this dipolar field are along this 
direction and X(k) is the dipolar interaction matrix el- 
ement for fc in the (110) direction (Fig. 3). In order to 
calculate the energy spectrum we now need to fix the 
size of the coherent dipole moment \fi\. According to our 
definition of the local mode the energy cost of flipping 
the direction of a single dipole out of the ordered ground 
state is defined to be Eq. This is equivalent to having 
2 \X(k = 0)| = Eo, which is the condition to have a gap- 
less mode at k — ► (§). Using this condition, the exper- 
imental value of Eo = 7KlKI determines the size of the 
coherent dipole moment as: ~ e • O.OlA. This value 
is indeed smaller than our previous estimation, which 
served as an upper bound on the size of the oscillating 
dipole moment (|2|). 

As we have proposed, the phase modulation in the 
(110) direction of the transverse atomic motion in the 
lattice, with energy E{k) (JsJ) should coincide with the 
Ti(110) phonon. In Fig. 4 we compare the experimen- 
tal values of Ti(110) taken from neutron scattering data 
with the calculated E(k), and we find that the agree- 
ment is excellent for all fc. From (||) and Fig. 3 we see 
that at the edge of the Brillouin zone the energy E(k) 
of the phonon should be just the bare energy of the lo- 
cal mode, Eq, since A(\/2vr/a) = 0. We also have that 
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at k = y/2n/a the dipoles have changed between the two 
configurations illustrated in Fig. 3, which are the two pos- 
sible configurations with alternating dipoles arranged on 
adjacent planes with the periodicity of the bcc unit cell. 

Since the empirical value of Eq that we used was taken 
from NMR data, the agreement we find with the phonon 
data taken from inelastic neutron scattering, emphasizes 
the self-consistency of our description. We stress that the 
value of Eq and the lattice vectors are the only empirical 
inputs used in the calculation, with the functional be- 
havior completely given by the lattice structure and the 
dipolar interactions. 

III. OFF-DIAGONAL-LONG-RANGE-ORDER 
AND CONDENSATION. 

We have found from the mean-field solution at zero 
temperature that the ground-state of the bcc phase con- 
tains a coherent-state of oscillating local-dipoles (Q). 
Since our method predicts the excitation spectrum of 
the Ti(llO) phonon with very good accuracy, we expect 
it to be valid at the finite temperatures for which the 
bcc phase exists. We therefore expect that the basic na- 
ture of the bcc phase will be well described by our re- 
sults, although the quantitative values may change due 
to the finite temperature. The coherent ground-state de- 
fines a global phase and breaks the gauge symmetry of a 
well-defined occupation number of local dipoles. In the 
limit k —> we find that the occupation number of the 
local-modes diverges as 1/fc, signaling macroscopic Bose- 
Einstein condensation in the zero-momentum state 

2 1 f E +X(k) \ lE /2 E / 

{Uk)=V {k) = 2 l^fcT " ) k -°2E(k) = ^k 

(9) 

where c is the sound velocity of the Ti(llO) phonon 
which is the natural excitation of the dipolar array. This 
is identical to the result for a Weakly Interacting Bose 
Gas (WIBG) problem solved by Bogoliubovfcj, where 
the divergence is related to the occupation of the zero- 
momentum state, i.e., the condensate fraction no/n 

( n k)wiBG = v ( fc ) ^fe-o — ^ (10) 

where in the WIBG case we have c as the k — > 
sound velocity, and = h 2 k 2 /2m is the free particle en- 
ergy. By comparing (|^) with (|l^) we find that in the bcc 
case the role of the condensate-fraction, the WIBG order- 
parameter, is taken by the parameter E , which is just 
2 |X(0)|. This can be seen directly from the form of the 
ground-state wavefunction (Q) where the pair-occupation 
is given by: 

=2u(*Mfc) = |^ (11) 



Equating the divergent part in ^ and (|o|) we can 
define an effective condensate fraction 

- = ^/2-^5 ±8% 

n mc z 10 

where we used for the velocity of sound c the values 
from our .calculation (||) (~ 130m/ sec) and from elastic 
constantsl!3 (~ 160m/ sec). It must be remembered that 
the mass m in (^|) is not necessarily the mass of a bare 
4 He atom since we are now dealing with condensation of 
local dipoles. Comparing with the condensate fraction 
at zero temperature in superfluid 4 HeL3, which is ~ 10%, 
we find that it is lower than the condensate fraction of 
the local modes in the bcc phase. We again mention 
that our result is for T=0 which can be depleted at finite 
temperature. 

It is clear that it is a non-zero coherent dipole moment 
fx that produces a dipolar interaction matrix X(k) which 
in turn implies finite pair occupation ( |TT| ) and a coherent 
ground-state. This is just the condition for the Bose- 
Einstein condensation of the dipoles in the bcc ground- 
state (|lj). We therefore have a broken gauge symmetry 
and a complex order-parameter in the form of the pair- 
occupation ( |Tl| ) . This function can be complex since the 
conditions on u(k) and v{k) allow for a relative complex 
phase between them, just as in the WIBG case. A similar 
condensation of local dipoles in all three orthogonal axes 
of local motion means that there are three independent 
phases at each lattice site, since orthogonal dipoles do not 
interact. The order-parameter in our case can therefore 
be described as a vector of three complex functions of 
independent magnitude and phase: 

/ |A**|e*»« \ 
$(r)= l^le^vW (13) 

If the cubic symmetry is not broken by external 
stresses, the magnitude of the coherent dipole moment 
in the three orthogonal directions should be the same: 
\Hx\ = \/j, y \ = \Hz\- I n the ground-state the phases are 
spatially uniform, while the excited state is described 
through a periodic phase oscillation, i.e. the Ti(110) 
phonon. The order parameter ( |l3| ) is to be contrasted 
with the order parameter of superfluid He, which also 
exhibits ODLRO and which has a single complex com- 
ponent. 

In the hep phase we do not expect the dipoles to or- 
der in a coherent state since the hexagonal geometry 
frustrates antiferroelectric-type configurations. Also the 
nearly isotropic potential of the hep lattice does not al- 
low the highly directional dipole moments as in the bcc 
case. Indeed there is good agreement between experi- 
ments and the-harmonic calculation of the phonons in 
the hep phaseEJ, indicating no strong quantum correc- 
tions, as in the bcc phase. 

The bcc 4 He is therefore a unique crystallographic 
phase having both Diagonal Long Range Order (DLRO) 



4 



of the solid lattice and Off-Diagonal Long Range Or- 
der (ppLBjQ4 of the local dipoles. It is not a 'super- 
solid 'EaOlla'tS in that it does not contain both a super- 
fluid and a solid, but is more similar to the superconduc- 
tors which have a DLRO of the atoms in-tlie lattice and 
ODLRO of the superconducting electronsEJ. This system 
is also distinct from the case of Bose-Einstein condensa- 
tion of a phonon mode which results in a static defomia- 
tion of the lattice and a structural phase-transitionEJ. 

Bose-Einstein condensation of local idefervts (vacancies) 
was previously considered for solid 4 Hdi3L3'lia, and is sim- 
ilar to our treatment. The main difference is that in our 
case the physical picture of the condensed local modes is 
not a local distortion of the lattice like a vacancy, and 
that the condensation i&jmiauc to the bcc phase. The 
estimate in these workst3'tlHl3 is that in the ground-state 
(T=0) the density of vacancies is ~ 10~ 4 per site (at mo- 
lar volume of 21cm 3 ), and a condensate density of ~ 10~ 7 
per site. In contrast we expect in the bcc phase a sizable 
fraction (10-30%) of condensed local-modes per site (12). 

The fact that the region of existence of the proposed 
supersolid phase in the phase diagram of solid 4 He should 
closely coincide with, the region occupied by the bcc 
phase, was shown inE3. In this work it was further shown 
that in the supersolid there should be a second-sound- like 
mode, which is an oscillation in the density of the local- 
defects. In our description of the bcc phase this suggests 
the possibility of an oscillation in the amplitude of the 
order-parameter <I>(r) (|l3|), that is in the amplitude of 
the coherently oscillating dipole moment. This is in con- 
trast to the Ti(110) phonon mode which is an oscillation 
in the phase of the order-parameter. This mode may be 
produced by modulating the density of the atoms so that 
the local excitation energies change and with them the 
amplitude of local motion and local electric dipole mo- 
ment. Unfortunately we do not expect such a mode to 
have measurable consequences which are different from 
the effects produced by usual longitudinal phonons. 

In concluding this section, we would like to mention the 
recent experiments on_the behaviour of implanted metalic 
ions (Cs) in solid 4 HeEH. These experiments are designed 
to look for evidence of time-reversal symmetry breaking 
which is equivalent to having a static electric dipole mo- 
ment. In our description of the bcc phase we do not find 
a static but a coherent-dynamic electric dipole moment. 
We point out that in these experiments a marked differ- 
ence between the hep and bcc phases has been found. 
In the bcc phase the electronic-spin relaxation of the Cs 
atom is extremely slow and this effect could be a result of 
the coherence and long-range order of the dipolar fields. 
The coherently oscillating 4 He electrons in the bcc phase 
will produce a very uniform electromagnetic interaction 
with the electronic spin in the Cs atom. By comparison, 
in the hep phase the spin polarization is extremely short 
lived, indicating a more random field environment. This 
result is in accord with our expectation that the coherent 
dipoles are unique to the bcc phase. Similar experiments 
in the future may allow a probe that will show directly 



the coherently oscillating-dipoles in the bcc ground-state. 

In these experimentsEll the hyperfine transition in the 
Cs atom was also measured. The energy shift of this tran- 
sition is sensitive to the shape of the confining cavity of 
the Cs atom inside 4 He lattice. The width of the transi- 
tion is a measure of the fluctuations in this cavity sizcE3, 
and the data show a much smaller spread in the bcc com- 
pared with the hep phase. Uncorrelated atomic motions 
of the 4 He atoms will increase the spread in instanta- 
neous cavity sizes due mainly to breating-like motion of 
the cavity walls (Fig. 5). This behavior is what we expect 
for the hep case. The correlated atomic motion in the 
bcc phase should result in a more constant cavity shape 
(Fig. 5) and a narrow signal, which is indeed measuredEil. 



IV. GROUND-STATE ENERGY AND THE 
STABILITY OF THE BCC PHASE 

The question of the relative stability of the different 
crystal structures in solid He has been a long standing 
one. The necessity for some non-Van-der-Waals interac- 
tions has been previously proposed to explain the occu- 
rancepQf fee over hep structure in the heavier rare-gas 
solidsEi The bcc phase is usually found to be more sta- 
ble than the close-packed hep-phase due to the large zero- 
point energy in the He solidsci The correlations between 
the dipoles in the ground-state that we have proposed, 
lowers the energy of the ground-state of the bcc phase and 
further stabilizes it with respect to the hep phase. The 
reduction in ground-state energy acheived by the coher- 
ent state of the dipoles-along one of the three orthogonal 
directions, is given byu3 



AE 



E 



E(k)-(E +X(k)) 
2E 



< 



(14) 



which is negative since X(k) < and E(k) < (Eq + 

Since the energy reduction integrand ( Km is non-zero 
only in the (110) directions it will give a small contribu- 
tion to the three dimensional phase-space integration. At 
zero temperature the summation in (14) will be confined 
to one dimensional sections along the (110) direction, so 
that the contribution will be zero. At the bcc temper- 
atures (~1.4K) the one dimensional chains in the (110) 
directions are broadened so that the summation in dlj) 
is now over finite volume sections of phase space. We can 
estimate the maximum width of the conical section in k- 
space as the momentum which corresponds to a Ti(110) 
phonon with energ y fc ^T, that is ~ 0.13A - 1 . The numer- 
ical integration of ( |14| ) over such volume sections gives an 
energy reduction of AE ~ — 2mK per atom. This result 
is in agreement with the experimentally interpolated en- 
ergy—difference between the bcc and hep phases of solid 
4 Heo, which is of the order of a few mK per atom. 

This reduction is less than 0.1 percent of the poten- 
tial and kinetic energies of the solid, and is therefore 
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very hard to calculate accurately theoreticallycJ. What 
is shown in the usual calculations is that the correlations 
between the motions of the atoms are essential in lower- 
ing the energy of the bcc phase, compared with the hep 
phase. Since part of the correlations in the atomic motion 
is described by our coherent dipolc model, we expect the 
condensation energy of the dipoles ( |l4| ) to be important 
in determining the stability of the bcc phase. At finite 
temperature the stabilization of the bcc phase compared 
to the hep phase is usually attributed to the lower zeca- 
point energy due to the lower T^llO) phonon energytJ. 
This is just the phonon which is softened by the long- 
range dipolar interactions that we have described, indi- 
cating again the importance of the coherent dipoles to 
the stabilization of the bcc phase of solid 4 He. Our pro- 
cedure may provide a good estimate of the small change 
in energy at the structural phase transition, by isolat- 
ing the degree-of-freedom which is most affected by the 
transition, i.e. the correlated atomic motion along the 
directions normal to the unit cell faces. 

The picture we propose is that the dipole condensation 
mechanism of the bcc phase competes with the lower po- 
tential energy of the hep phase due to its higher coordina- 
tion number. If the hep phase has a large enough volume 
(through thermal expansion or introduction of 3 He im- 
purities) its potential energy is increased until a critical 
point is reached where the bcc phase has a lower total 
energy due to the dipolar-condensation-.energy reduction 
(|l4|), which is absent in the hep phasea. At this critical 
point the structural phase transition occurs. 

By comparison, solid 3 He has a stable bcc phase due to 
the larger kinetic energy of this lighter isotope. This in- 
creased zero-point energy causes the less dense bcc phase 
to have a lower ground-state energy than the hep phase 
even at T=0. Since the 3 He atoms are fermions with a 
spin 1/2 nucleus, the oscillating electric dipoles are not in 
resonance and can not be treated as bosonscl. We there- 
fore do not expect a coherent state of the atomic mo- 
tion as in the bcc 4 He, and bcc 3 He is stable due to its 
large zero-point kinetic energy alone. On the other hand, 
at low enough temperatures where the bcc 3 He becomes 
an antifcrromagnet, there could be correlations involv- 
ing both the nuclear spin and electric dipole degrees of 
freedom. 



V. FERMIONIC EXCITATIONS 

In addition to the fluctuations of the phase of the co- 
herent dipole ground-state (i.e. Ti(llO) phonons), there 
can be a localized 'flip' of a dipole so that it is in anti- 
phase (phase difference of ir) relative to the rest of the 
dipoles, in the ground-state configuration. Such an ex- 
citation is naturally treated as a Fermion since such a 
flipped dipole is antisymmetric with respect to the other 
dipoles, that is with respect to the global phase 8 ( |l3| ) in 
one (or more) of the orthogonal directions of local motion 



[x, y, z). The flipped dipole is no longer a dipolar image 
of the nearest-neighboring dipoles but an anti-image, and 
will be treated with Fermi-Dirac statistics. 

An anti-phase localized-mode (a fermion) is not part 
of the correlated ground-state, but nevertheless will feel 
the effect of the Bose excitations (Ti(llO) phonons) of 
the dipolar array as they interact with it. The effective 
Hamiltonian describing such a fermion should therefore 
contain a term describing the creation and anihilation 
of pairs of fermions from the ground-state by a phonon 
(Boson). This is an off-diagonal term that describes the 
fluctuation caused by a Ti(llO) phonon of energy E(k): 
it changes a fermion 'particle' into a 'hole' and vice versa. 
The terms 'particle' and 'hole' are with respect to the 
ground-state which has occupation of pairs of localized- 
modes (i.e. not an 'empty' vacuum). 

In addition there should be a term that describes the 
excitation energy of the bare fermionic localized-mode, 
that is Eq. This is just the energy to 'flip' a dipole 
from the ground-state so that all it's interactions with 
the other dipoles of the ground-state change sign, i.e. 
-2X(0) = E . 

The many-body effective Hamiltonian that we there- 
fore propose is 

H D = E(k) (c\cl k + c fe c_ fe ) - K (4c fc c! fc c_fc) 

k k 

(15) 

where ct , c k are the creation and annihilation operators 
of the anti-phase (Fermionic) localized-mode. The first 
term in ( |l5| ) is the 'kinetic' term due to the phonon-roton 
branch, where the localized-modes are created/anihilated 
in pairs. The energy E(k) is the energy of the Ti(llO) 
phonon excitation (0). In addition there is a finite 'po- 
tential' energy if there is a finite density of unpaired 
fermions, which is Eq. In the absence of the second term 
we have just the Bose ground-state written in terms of 
fermionic pairs. 

We linearize the equations of potion that follow from 
(H), similar to the BCS method^ 

i% c k = -E(k)J_ k + A fe c fe i% c_ k = -E(k)c k - A£cl & 

(16) 

where we used the Fermi anti-commutation relations: 
{cfc,c|, j = 1 { c *:i c -/c} — 0, and we define 

A k = A* k = E = J2 v k (cW) , ( c I c fe) = hV k = E 

k k 

(17) 

From (|l7|) we see that the symbol Eq will now indicate 
a finite density of fermions. 

The equations of motion have the following eigenval- 
ues: 



G 



E f (k) - E Q -E(k) 
-E(k) E f (k)+E 



Q^E f {k) = x lEl+E{kf 
(18) 



We can now solve the equations using the Rogoliubov- 
Valatin transformation for superconductivity^: 



c k = u k a k + Vka_ k cl fe = ~v k a k + u k a f _ k 



with the functions u k ,v k given by 



4 = ~H 



E 



Ef(k)," ' 
The ground state is 



E 
Ef(k) 



(19) 



(20) 



a k |0) = => |0) = [Ja-fcQt \vac] 

k 

= n [u k + v k c\J_^j \vac) 

k 



(21) 



In Fig. 6 we plot the energy spectrum ( |l8| ) compared 
with the other phonon modes in the (110) directior£3. It 
is clear that this optic-like branch should be detectable 
in the low momentum range where it is not masked by 
the signal from the accoustic phonon modes. There is 
at present no high resolution neutron-scattering data in 
this energy and momentum range, and this prediction 
can be hopefully checked in future experiments. This 
mode could also be observed by Raman scattering, as a 
peak at energy Eq. 



VI. SCATTERING INTENSITY 

A neutron scattering inelastically from the solid will 
create/anihilate an elementary excitation. An excitation 
from the effective ground-states of the Ti(110) phonon 
(0) and of the fermionic mode ( ^l| ) involves an anihilation 
of a pair of local-modes, leaving an unpaired local-mode. 
We therefore expect the experimcntaly measured neutron 
scattering intensity to be proportional to the density of 
local-mode pair-occupation at each wavevector k. This 
gives us the following results for the two modes: 

Ti(110) phonon: 



E n 



E{k) 



E(k)V 
E J 



Fermionic excitation: 

1 E{k) 



I oc 



2E f (k) 



1 



(E /E(k)j 2 + 1 



(22) 



(23) 



Both functions fl22p,(23) are plotted with arbitrary 
scale in Fig. 7. We sec that the intensity of the Ti(110) 
phonon is such that at small k it behaves as 1/k 



which is typical for phonons at low k and was seen 
experimentalljO, but goes identically to zero at the edge 
of the Brillion zone where E(k) — > Eq. The fermionic ex- 
citation has an opposite behavior by increasing linearly in 
intensity with k, until it saturates at the edge of the Bril- 
lion zone. The expression for the intensity of the Ti(110) 
phonon is similar to the expression of the intensity of the 
phonon-roton excitation spectrum of superfluid 4 HeEZI, 
where it agrees very well with the experimental results. 

These predictions for the bec phase have yet to be 
checked experimentally. 



VII. CONCLUSION 

In this work we have investigated the nature of the 
quantum correlations in the bec phase of solid 4 He. We 
identified a three component complex order parameter 
and Bosc-Einstein condensation in this phase, though not 
a 'super-solid'E3, i.e. no superfluid component. There can 
be further manifestations of the ODLRO of the dipolcs 
in the bec phase which we have not explored yet, such 
as macroscopic topological defects in the complex order- 
parameter. The order-parameter or condensate-fraction 
can also serve as an extra thermodynamic variable, and 
this opens the possibilty of more complicated internal 
dynamics in the bec solid, such as the phenomenon of 
second sound in superfluid 4 He. 

We predict that a local excitation of a dipole out of the 
coherent ground-state will behave as a Fermion, and we 
calculate its energy spectrum. We find it to behave as an 
optical-like branch in the (110) direction. Finally we cal- 
culate the scattering intensity as a function of wavevec- 
tor k for both the Bose (Ti(110) phonon) and Fermi 
(new optical-like branch) excitations. All these predic- 
tions await high-resolution neutron and Raman scatter- 
ing experiments to be compared with. 
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VIII. APPENDIX A: COMPARISON OF BOSE 
EXCITATIONS WITH KLEIN-GORDON 
HAMILTONIAN. 

We would like to point out that the Hamiltonian de- 
scribing the localized dipoles (|5|) is similar to the Klein- 
Gordon (KG) Hamiltonian foj; a single spinless boson, 
written in its first-order formc3: 

Klein-Gordon: 

H K g = £k (<t z + iv y ) + mc 2 a z + e$T (24) 
Localized dipoles (||): 

Hi oc = X(k) {a z + ioy) + E cr z (25) 

where Uj are the Pauli matrices, £k = p 2 /2m, e is the 
electric charge, $ is the electrostatic potential, m is the 
KG-particle's mass and c is the velocity of light. 

We have written the dipolar hamiltonian ( p5[ ) in the 
basis of a two component wavefunction 



becomes the equation for X(k) in the dipolar case: 



C-k 



(26) 



In this representation we see that exciting a local dipole 
out of the ground-state configuration (c k ) has bare energy 
Eq while destroying an excited dipole has minus this en- 
ergy. There is a freedom of choice weather to define the 
positive excitation to be a flipping of an up dipole to a 
down dipole or vice versa. The sign of the energy of the 
dipolar bosons therefore represents this freedom which 
corresponds to two equivalent dipolar configurations with 
a 7r phase difference. 

The two-component wavefunction of the KG Hamilto- 



- U + — V>° 
2 \ mcr 



mcf 



(27) 



where ip is the original wavefunction of the second- 
order KG equation, and V>° = (-§- t + f $) ij>. The Hamil- 



tonians (|24j,[25|) are the similar except that the KG density 
is not normalized to 1 but to (p) = E/mc 2 , describing 
the relativistic increase in the density with velocity. By 
comparing the two Hamiltonians ( |25| , p4] ) we identify that: 
Eq «-> mc 2 , X(k) <-> Efc, which gives the equivalence of 
the two hamiltonians. 

The peculiarities of the KG equation appear when 
there is a potential V — e$ (the Klein paradox for exam- 
ple). The equation for the momentum of the KG prticle: 



2mc z e k = h 2 c 2 k 2 = (E(k) - Vf 



^k = 



\]{E(k) - V) 2 - {mc 2 

he 



(28) 



X(k) 



(E(k)-V) 2 -(E(k) ) 2 

2En 



(29) 



We see from (|29|) that there is a region of energies 
where the interaction parameter X{k) is positive and 
a region where it is negative. We saw above that the 
condensation of the dipoles in the bec phase is charac- 
terized by a negative X(k) which also gives a gapless 
excitation spectrum at k — > 0. The excitations with 
E(k) > V + Eo,E(k) < V - E are therefore not con- 
tributing to the coherent long-range order. A fermionic 
excitation is a local destruction of the coherent order , 
and indeed costs at least Eq (for the free case with V = 0) 
to create (|T8|). 

In the case of the KG equation the sign of the enrgy 
indicates the charge of the particle/ antiparticle, which 
have oposite charges. Charge conjugation therefore in- 
terchanges between the two. What is the meaning of the 
different signs of the energy of the dipolar bose excita- 
tions in the bec case ? From our definition of the second- 
quantized description of the dipoles, the meaning of the 
sign of the energy is that the field of resonating localized- 
dipoles can have two global configurations shifted by ir 
(Fig. 2). These two configurations are identical with re- 
spect to the magnitude of the energy spectrum, but in 
each the operation of spin flip changes from up^down to 
down— > up. We can therefore identify two "charges" for 
the bec to distinguish between the two shifted phases. 
Further we find that as in the KG case the operation 
of charge conjugation (which reverses the signs of the 
dipoles) moves us between the two solutions. 



IX. APPENDIX B: SYMMETRY BREAKING OF 
THE FERMIONIC EXCITATIONS. 



We see from ( 20 ) that when there is no fermion present 
(i.e. if we put Eq = in (|L7|) ) the ground state has equal 
numbers of fermions and holes. The symmetry between 
particles and holes is broken by the free fermion quasi- 
particle (or quasihole), and the sign of the symmetry- 
breaking parameter Eq determines which of the two kinds 
is present. The hole/particle are with respect to the equi- 
librium occupation by pairs of fermions in the ground- 
state. 

A single flipped dipole described as the fermionic ex- 
citation, breaks the symmetry between the number of 
up/down dipoles and creates a residual globally oscillat- 
ing dipole moment. The parity P symmetry with respect 
to reflection along the axis of the global dipole (let us 
choose to be z) is broken. The charge C symmetry is also 
broken since the direction of the global dipole is flipped 
under charge conjugation. The time reversal symmetry 
T is unbroken since the oscillating globel dipole does not 
define a unique time direction. We therefore have that 
the global CPT symmetry is preserved, as is the CP and 
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T symmetries individually. The symmetry-breaking pa- 
rameter in (|l5| ) is the sign given to Eq, which corresponds 
to choosing an up or down dipole to flip. In second quan- 
tization langauge this is the choice between an unpaired 
particle or hole. The broken symmetry is not of the U(l) 
group such as the $ ( |l3| ) order parameter, but has a Z(l) 
discrete symmetry. 

We now compare this with the situation of the two- 
dimensional massive Dirac particleO. The Hamiltonian 
describing a single fermionic excitation ( |l5| ) can be writ- 
ten as 



at V c ] _ k 



E(k) 

Eq 



1 

1 

1 

-1 



c -k 

Ck 
,t 

'-fe 



while the 2D Dirac particle is described by 



in— 



dt\x 



2^ / <P 

+ mc a z 

X 



(30) 



(31) 



By assuming momentum p in the x-direction only we 
write 



dt\x V -ibdx I \ x 



mc 



1 \f<p 



o -l A x 



(32) 



where a — (a x ,a y ,d z ), at the Pauli matrices, and tp, x 
are the particle/antiparticle scalar wavefunctions. There 
is now complete analogy between (|3C]) and (|2|). The 
symmetry-breaking parameter Eq is identical to the mc 2 
parameter in the 2D Dirac equations-of- motion (^) . The 
symmetry that is broken by choosing a non-zero mc 2 Jui 
two-dimensions is the time-reversal symmetry (TRS)t3. 
As shown in ( |3l| ) , the time component of the momentum- 
energy vector in two dimensions is taken by the z axis. 
The broken parity P in the z-axis for the fermionic ex- 
citation of the bcc phase is here replaced by the TRS 
breaking of the heavy two-jiimensional Dirac particle. 

Similar to the AndersonE3 transformation of the BCS 
problem to a magnetic Hamiltonian we can transform 
© using: 

n k = c\c k 4 c -fc = a k / 2 c fcC-fc = (Tfe /2 

=> n k ri-k = - (of + 1) 4 c -fc + c k c -k = of (33) 



This Hamiltonian describes a fictitious magnetic field 
acting on the spin 



B 



1 



s k x - 



-E z 



(35) 



where we replaced the potential energy with the con- 
stant V k = E @. The magnetic field @ can be 
compared with the BCS resultE^I 



(36) 



The alignment of the spins in the ground-state is shown 
for the two Hamiltonians in Fig. 8. In the BCS problem 
the sign of the symmetry-breaking field in V (|3^) has to 
be positive so that it induces ferro-magnetic interaction 
between the fictitious spins, otherwise there will not be 
any rotation of the spins across the Fermi-energy. Only 
when the spins rotate do they go through the point where 
the spin is entirely in the iry-plane. At this point the 
state has no defined occupation number but a well defined 
phase, while on both sides of the domain-wall there is well 
defined occupation and no phase. This superconducting- 
phase at the fermi energy is just the angle of the fictitious 
spin in the xy-plane, and there is broken U(l) symmetry 

In the Dirac case the symmetry-breaking field E in 
( |35| ) is a constant external field. It can have both signs, 
which control the direction along the z-axis that the ro- 
tated spin has in the middle of the domain-wall. This 
spin describes weather a particle or a hole is occupied, 
while away from the k = point the spins are in the 
xy-plane, with a well defined phase. The symmetry that 
is broken is therefore the binary Z(l) (±) symmetry 

In comparison with the BCS problem we see that in the 
fermionic excitation in the bcc the symmetry-breaking 
parameter is a finite density of unpaired fermions: 

(c\ck*J 7^ (|l7|). The ground-state without unpaired 

fermions is a 'vacuum' of pairs of particle-holes in equal 
numbers. In the BCS problem the symmetry-breaking 

parameter is a finite pair-density: ( c \^-k) ^ ^- The 
ground-state in the absence of electron pairing is just a 
finite density of electrons below the fermi-energy and zero 
above. In this respect the two problems are 'complemen- 
tary'. 



where the crjj. are Pauli spin-1/2 operators. The basis 
is such that an up-spin in the z-direction represents an 
empty pair, while a down-spin represents an occupied 
pair. 

The resulting Hamiltonian is: 



Hmaa — 



5l> 



1) 



(34) 
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FIG. 1. The potential- well of an atom in bcc 4 He along dif- 
ferent directions. The energy difference Eo between the low- 
est two energy levels (dashed lines) are: (111)- 59. 5K, (110)- 
27.6K, (100)- 10.6K. 




FIG. 2. The coherent dipole arrangement in the 
ground-state of the bcc phase, oscillating between these two 
configurations. Dipoles with same phase have same shade. 
The sum of the dipole-dipole interaction (Eq.3) for a unit 
dipole-moment is: -0.08 (A 3 ). 
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FIG. 3. The calculated interaction matrix X(k) (Eq.4) as 
a function of the wavevector k, for the three phonon modes 
that could affect a dipolar array. The dipole moment has been 
normalized to give a gapless mode: X(k = 0) = —Eo/2. The 
unit-cell dimensions (A): a = 4.12/2, a' = a\/2. Also shown 
are the two arrangements of the dipoles in the extreme points 
along the oscillation in the (110) direction. 
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FIG. 4. The experimental data [13] (solid circles) for the 
Ti(110) phonon compared with the calculation (Eq.8) (solid 
line). Also shown is the energy of the bare local-mode 
E = 7K. 
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FIG. 5. Schematic cavity shapes for a Cs atom (grey circle) 
inside 4 He solid (empty circles with arrows), (a) Correlated 
atomic motion in bec: constant cavity shape, (b) uncorrelated 
atomic motion in hep: randomly fluctuating cavity shape. 



11 




Wavevector K 



* 



0.2 



0.8 



1 



0.4 0.6 
Wavevector (A -1 ) 

FIG. 6. The spectrum of the fermionic optic-like mode 
(Eq.18, solid line) compared with the experimentally mea- 
sured phonons in the (110) direction [13]. 
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FIG. 7. The calculated intensity of the inelastic neutron 
scattering by the Ti(110) phonon (Eq.22, solid line) and by 
the new fermionic mode (Eq.23, dashed line). 





FIG. 8. Arrangement of the spins in the mag- 
netic-analogue-Hamiltonian. BCS: (a) without pairing, (b) 
with pairing. Dirac: (c) without unpaired fermion, (d) with 
unpaired fermion. 
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